The cascades route to chaos 
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The presence of a period-doubling cascade in dynamical systems that depend on a parameter is 
one of the basic routes to chaos. It is rarely mentioned that there are virtually always infinitely many 
cascades whenever there is one. We report that for one- and two-dimensional phase space, in the 
transition from no chaos to chaos - as a parameter is varied - there must be infinitely many cascades 
under some mild hypotheses. Our meaning of chaos includes the case of chaotic sets which are not 
attractors. Numerical studies indicate that this result applies to the forced-damped pendulum and 
the forced Duffing equations, viewing the solutions once each period of the forcing. We further 
show that in many cases cascades appear in pairs connected (in joint parameter-state space) by 
an unstable periodic orbit. Paired cascades can be destroyed or created by perturbations, whereas 
unpaired cascades are conserved under even significant perturbations. 



In Fig. [T] as /i increases towards a value /i^? 3.57, one 
encounters a family of periodic orbits that undergo an 
infinite sequence of period doublings with the period of 
these orbits tending to oo. The appearance of infinitely- 
many such (period-doubling) cascades is one of the most 
prominent features observed in the study of parametrized 
maps. Cascades were first reported by Myrberg in 
1962 [1] (c/. Fig. [T]), and Robert May popularized their 
existence to a huge scientific audience [2j . They are found 
in a large variety of contexts: in Raleigh-Bernard convec- 
tion [3], damped bouncing balls [4], Van der Pol oscilla- 
tors (c/. Fig.[2|, vibratory ball milling [F, dust charge 
fiuct nation in plasma systems [7^, external optical injec- 
tion in lasers [8] , delay oscillators [12] , vibrating damaged 
structures [TT, glow discharge [TT, combustion [15 , very 
slow classical Cepheid stars [16 , neuron and pancreatic 
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FIG. 1: Cascades and their connected components. 

The attracting set for the logistic map F{/j.,x) = /ix(l — x) is 
shown in blue. There are infinitely-many cascades, each with 
infinitely- many period-doubling bifurcations. Each saddle- 
node bifurcation creates both a cascade and a path of un- 
stable orbits, the latter shown as black curves for periods up 
to six. We call the paths of unstable orbits the stems of the 
cascades. Below the largest gaps or windows in the blue re- 
gions, we give the period of the associated cascade. For this 
particular map, the black curves of unstable orbits extend to 
/J = +(X). Here, at each sufficiently large /x, there is exactly 
one orbit {k points) in each period-Zc stem. 




FIG. 2: Cascades in the double-well DufRng equation. 

The attracting sets (in black) and periodic orbits up to period 
ten (in red) for the time-27r map of the double-well Duffing 
equation: x'\t) + 0.3x\t) - x{t) + {x{t)f + {x{t)f = /xsint. 
Numerical studies show regions of chaos interspersed with re- 
gions without chaos, indicating that our Off-On-Off Chaos 
Theorem applies to this situation. 



cells [17] , dripping faucets [T8^ , bouncing droplets on soap 
films p!9l, Belousov-Zhabotinsky reactions [20^, and in 
the cromorn, a medieval musical instrument [2T] . 

We find that for many systems depending on a pa- 
rameter, every periodic orbit is part of a cascade; for 
such systems, cascades are as fundamental as periodic or- 
bits themselves. The scaling properties of individual cas- 
cades have been studied for cascades in nearly quadratic 
maps [22], but there are only a few results about the 
existence of cascades [23]. Furthermore, the mathemati- 
cal and scientific literature focuses on the study of single 
cascades. In this Letter, we describe the results of a new 
general theory of cascades, which explains why cascades 
exist and why chaotic dynamical systems often have in- 
finitely many cascades. 

Two kinds of periodic orbits meet at period- 
doubling bifurcations. For a map F(/i, x) that 
depends on a parameter /i, a point {fi^xo) is a 
period-p point if F^(/i,xo) = xq, where p > 
is chosen as small as possible. Writing x^+i = 
F{fi^Xn), its (periodic) orbit is the set of points 
{(/i, xo), (/i, xi), (/i, X2), • • • , (/i, Xp-i)}. By the eigen- 
values of that orbit, we mean the eigenvalues of the 
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FIG. 3: Bounded paired cascades. The attracting set 
for the function F(/j,x) = h{/j)x{l — x), where = 
/i(1.18 + 0.17cos(2.4/i)). Cascades Ci and C2 constitute a 
pair of cascades, since they are connected by a family of or- 
bits (black curves). Cascade C3 is unbounded and unpaired, 
as indicated by the black curve continuing to the right to 
/i = 00. 



Jacobian matrix D^F^^ii^xq). If the map is one- 
dimensional, then by the derivative of the orbit, we mean 
the derivative of F^. A bifurcation orbit is any orbit 
having an eigenvalue with absolute value 1. 

We call a periodic orbit a flip orbit if the orbit has an 
odd number of eigenvalues less than -1, and -1 is not an 
eigenvalue. (In one dimension, this condition is: deriva- 
tive < — 1. In dimension two, flip orbits are those with 
exactly one eigenvalue < — 1.) All other periodic orbits 
are called regular. A family switches between flip and 
regular orbits at a period-doubling bifurcation orbit since 
this bifurcation orbit always has an eigenvalue equal to 
-1. 

For il) in some interval (a, 6), let be a path 

of re^^/Zar periodic points in (/i, y) space and write [1^('0)] 
to denote the periodic orbit of the periodic point Y{tlj). 
We say (/i('^), [F(V^)]) is a cascade if the path contains 
infinitely many period-doubling bifurcations and has or- 
bits with all the periods {/c, 2/c, 4/c, 8/c, • • • } for some pos- 
itive integer k. 

Assumptions. All results in this paper assume that 
the maps y 1-^ y) are infinitely differentiable jointly 
in the parameter ja and in phase space y and that all 
of the periodic orbit bifurcations are generic^ meaning 
that every bifurcation orbit is a standard nondegenerate 
period- doubling, saddle-node, or Hopf bifurcation. For 
example, generic maps have neither symmetries nor sym- 
metry breaking bifurcations. Nor do they have pitchfork 
or period-tripling bifurcations. We also assume that for 
any parameter interval [/ii,/i2], all periodic orbits are 
contained in a bounded set in phase space. 

The cascades route to chaos in dimensions one 
and two. Usually the term "route to chaos" refers to 
formation of chaotic attractors, in which case there are 
several possible routes to chaos. In this Letter, we con- 
sider routes from a parameter having no chaos to one 
with chaotic sets^ where the sets are not necessarily at- 
tractors. In fact, we require only one aspect of chaos: we 



say that there is chaos at a particular /i if there are in- 
finitely many regular periodic orbits. For example, there 
is chaos whenever there is a transverse homoclinic point. 
That is equivalent to having a horseshoe for some iter- 
ate of the map, and it implies there are infinitely many 
regular saddles in two dimensions, and infinitely many 
regular unstable orbits in dimension one. This definition 
of chaos is sufficiently general as to include having one 
or multiple coexisting chaotic attractors, as well as the 
case of transient chaos. We say a map has no chaos at 
a particular ja if there are at most finitely many regular 
periodic orbits. 

Cascades Route to Chaos Theorem. Assume the 
parametrized map F{ii^ •) has a one- or two-dimensional 
phase space. If at parameter value jii there is no chaos, 
and at 112, there is chaos and there are at most finitely 
many attracting orbits ( and finitely many repelling orbits 
in dimension 2), then F has infinitely many cascades be- 
tween jii and /i2- 

Thus for one- and two- dimensional smooth families, 
the only route to chaos is through infinitely many cas- 
cades. It is for instance impossible to get chaos at 112 
from a single cascade. 

If we were to omit the assumption that there are in- 
finitely many regular orbits, the conclusion would be 
false. For example, in the quadratic map let fip denote 
the Feigenbaum parameter value, i.e., the first parame- 
ter value where the period doublings accumulate. At jhf 
there are infinitely many periodic orbits, and all but one 
are flip orbits. Furthermore /ii? is preceded by a single 
cascade. 

Each cascade will have some minimum period p and 
will have periodic orbits of periods p, 2p, 4p, 8p, • • • . If 
the map is one-dimensional or is two-dimensional and 
dissipative in the sense that there are no repelling peri- 
odic orbits, then each of these periods is the period of 
some attracting orbits in the path, though not all orbits 
will necessarily be attracting. 

Based on numerical studies, a number of maps appear 
to satisfy the conditions of the Cascades Route to Chaos 
Theorem. Note that these numerical verifications involve 
significantly more work than just plotting the attracting 
sets for each parameter, since we are concerned about 
both the stable and the unstable behavior to determine 
whether there is chaos. Examples include the time-27r 
maps for the double- well Duffing (Fig. |2|, the triple- well 
Duffing, and forced-damped pendulum (Fig. [4|, as well 
as the Ikeda map used to describe the field of a laser 
cavity, and the pulsed damped rotor map. 

We now discuss the creation of a path (/i('0), [F(V^)]) 
of regular orbits. Consider the procedure of starting at 
a regular periodic orbit and following the path of peri- 
odic orbits containing it in (/i, ?/)-space, following only 
regular orbits. Follow the path - such as numerically - 
through saddle-node bifurcations by reversing the direc- 
tion in the parameter space. When period-p orbits reach 
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FIG. 4: The forced-damped pendulum. For this figure, 
periodic points with periods < 10 were plotted in red for 
the time-27r map of the forced-damped pendulum equation: 
x" (i) + 0.2x'(t) + sin(x(t)) = a;cos(t), indicating the general 
areas with chaotic dynamics for this map. Then the attracting 
sets were plotted in black, hiding some periodic points. Much 
more detailed calculations confirm that interspersed with the 
chaos, there are some parameter ranges without chaos. 

a period-doubling bifurcation (/io, \X\)^ follow the period- 
2p orbits. They are regular orbits, whereas the extension 
past (/io, \X\) of period-p orbits are flip orbits. Likewise, 
at a period-halving bifurcation, follow along the branch 
of regular orbits. Following a path away from a cascade 
can lead to ^ oo or it can lead to a second cascade. 
We call an entire connected path of regular orbits a com- 
ponent. Two orbits are in the same component if and 
only if it is possible to follow a path of regular orbits 
from one orbit to the other continuously, without jumps, 
though the period is allowed to double many times. 

There are two types of cascades. Only one type 
of cascades occurs for quadratic maps, but there is a sec- 
ond type of cascade as well. Namely, the component of a 
cascade either lies in a bounded region of (/i, ?/)-space, or 
it does not. Accordingly we call the cascade bounded 
or unbounded, respectively, though the adjective re- 
ally describes the cascade's component. Each cascade in 
the quadratic map (c/. Fig. [T]) is unbounded because a 
branch of periodic orbits extends from each cascade to 
/i = oo. The constant period branch that extends to 
infinity is called the stem of the (unbounded) cascade. 
Thus the theory of quadratic maps is the theory of un- 
bounded cascades. We call the period of the orbits in the 
stem, the stem period of the unbounded cascade. 
(We only treat the case of an unbounded cascade which 
has a constant period stem.) 

We refer to two cascades as a pair if they are in the 
same component; that is, if the two cascades have a path 
of periodic orbits running from one the other (c/. Fig.[3|. 
We have shown the following rigorous result. 

Bounded Paired Cascades Theorem. Bounded 
cascades always come in pairs. 

The theorem below shows that paired cascades are 
common. Once again we consider x in a one- or two- 
dimensional phase space. We say that map F has off- 
on-off chaos for /ii < /i2 < Ms if there is no chaos at fii 
and /i3, whereas at /i2 5 F is as assumed in the Cascades 
Route to Chaos Theorem. 

Off-On-Off Chaos Theorem. // F is two- 
dimensional and has off- on- off chaos for /ii < /i2 < MS; 



then F has infinitely many bounded paired cascades and 
at most finitely many unbounded cascades. 

Our numerical studies indicate that the time-27r maps 
of the forced double-well Duffing (Fig. |2| and forced 
damped pendulum (Fig.[4| have off-on-off chaos, and that 
it occurs on multiple non-overlapping parameter regions. 
Here 27r is the forcing period. In fact we cannot prove 
there are only finitely many periodic attractors, though 
we find very few. These systems have no periodic re- 
pellers. 

Conservation of unbounded cascades. If in addi- 
tion to the hypotheses of the Cascades Route to Chaos 
Theorem, we assume further that {i) at /ii there are no 
regular periodic orbits and {ii) at /i2 there are no attrac- 
tors (nor repellers in dimension two) , then the component 
of each regular orbit at /i2 contains a cascade that is be- 
tween /ii and /i2, and distinct regular orbits at /i2 are in 
different components. Starting at each regular orbit at 
/i2, there is a unique path of regular orbits - initially in 
the direction of jui - ending in a cascade. Thus the num- 
ber of cascades hitting the boundary of the parameter 
interval depends only on the number of regular orbits at 
/i2, independent of the behavior of the map between fii 
and /i2- 

This leads to a heuristic conservation principle for 
unbounded cascades: if F{jj.^x) satisfies property {i) 
for /i sufficiently negative and property (ii) for /i suffi- 
ciently positive, then any perturbation F{ii^x) + g{/J^^x) 
such that F dominates ^ for ^ oo has the same set 
of unbounded cascades. 

Three rigorous examples of this idea are encapsulated 
in the following one-dimensional maps 

F(/i, x) = /i — + ^(/i, x) (quadratic), (1) 
F(/i, x) = jj.x — x^ -\- g{jj.^ x) (cubic), 
F(/i, x) = x^ — + + g{fi, x) (quartic), 

where for some real positive 

|^(/i,0)| < P for all /i, and (2) 
\gx{/i,x)\ < P for all 

For such ^, each of the three maps each has no regular 
periodic orbits for ja sufficiently negative, and for ja suffi- 
ciently large has a one-dimensional horseshoe map. The 
quartic map was chosen so that when g is identically 
and /i > 0, the graph has two minima, at (±y7^, — /i^/2), 
and the local maximum is at (0, +/i^/2). The conditions 
on g guarantee that it does not significantly affect the 
periodic orbits when is sufficiently large; in particular 
it does not affect their eigenvalues, so it does not affect 
the number of period-p regular periodic orbits for large 
Hence one can check that all three maps have no 
regular periodic orbits for /i very negative and have no 
attracting periodic orbits for fi very positive, and for suf- 
ficiently large all periodic orbits are contained in the 
set [—2y/JI,2y^]. Thus the following holds. 
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56 


330 
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156 


1152 
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14 


400 


4064 


9 


28 


1092 


14560 


/O 10 


- 2^ /2k 


- 3^ /2k 


- 4^ /2k 



TABLE I: The number of unbounded cascades of stem period 
k for any large-scale perturbation of the quadratic, cubic, and 
quartic maps, respectively labeled Quad(k), Cubic(k), and 
Quart{k). 



Conservation of cascades theorem. For any F of 

the functions in Eqn. [7j with g chosen as in Eqn. [H the 
number of stem-period-k unbounded cascades is indepen- 
dent of the choice of g. 

The number of unbounded cascades in each of these 
three cases is summarized in our Table. The relative 
sizes reflect the complexity of the three horseshoe maps; 
the quadratic map has the complexity of a standard two- 
branch horseshoe, whereas the cubic has a horseshoe with 
three branches, and the quartic has four branches. 

The conservation principle works for two- and higher- 
dimensional maps as well. We have shown in [25 that 
large-scale perturbations of the two-dimensional Henon 
map conserve unbounded cascades. In [24 , we have 
shown that there is conservation of unbounded cascades 
for a coupled system of TV quadratic maps. 
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